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Abstract: We analyse the interacting theory of charged fermions, scalars, pseuso- 
scalars and photons propagating in 2-dimensional curved spacetime in detail. For 
certain values of the coupling constants the theory reduces to the gauged Thirring 
model and for others the Schwinger model incurved spacetime. It is shown that 
the interaction of the fermions with the pseudo-scalars shields the electromagnetic 
interaction, and that the non-minimal coupling of the scalars to the gravitational 
field amplifies the Hawking radiation. We solve the finite temperature and density 
model by using functional techniques and in particular derive the exact equation of 
state. The explicit temperature and curvature dependence of the chiral condensate 
is found. When the electromagnetic field is switched off the model reduces to a 
conformal field theory. We determine the physically relevant expectation values 
and conformal weights of the fundamental fields in the theory. 



1. Introduction 

The response of physical systems to a change of external conditions is of em- 
minent importance in physics. In particular the dependence of expectation 
values on temperature, the particle density, the space region, the imposed 
boundary conditions or external fields has been widely studied [1]. Despite 



2 A. Dettki, I. Sachs and A. Wipf 

all these efforts we are still unable to understand, for example, the mecha- 
nism leading to the spontaneous symmetry breaking of the SUAiN) in low 
temperature QCD [2]. Clearly such subtle effects require a better under- 
standing of the nonperturbative effects and in particular nonperturbative 
vacuum sector of gauge theories. From our experience with 2-dimensional 
gauge theories [3] which we suppose to mimic one-flavour QCD [4], we are 
lead to believe that gauge fields with windings are responsible for the non- 
vanishing chiral condensate and in particular its temperature dependence. 
A related problem is how quantum systems behave in a hot and dense en- 
vironment as it exists or existed in heavy ion collision, neutron stars or the 
early epochs of the universe [2] . 

On another front there has been much effort to quantize selfinteract- 
ing field theories in a background gravitational field [5]. For example, one 
is interested whether a black hole still emits thermal radiation when self- 
interaction is included. Due to general arguments by Gibbons and Perry [6] 
this question is intimately connected with universality of the second law of 
thermodynamics . 

Rather than seeking new partial results for more general and realistis- 
tic 4-dimensional systems we have chosen an idealized 2-dimensional model 
with self-interaction to investigate the questions mentioned and others. It is 
a theory containing photons-^, charged massless fermions, scalars and pseu- 
doscalars in interaction with themselves and a gravitational background 
field. The model has the action 

where F^i^ is the electromagnetic field strength, the gamma-matrices in 
curved space are related to the flat ones as 7'^ = e'^7'*, V ^ = d^+iuj^—ieA^ 



^ Altlioutli pliotons in 1-1-1 dimensions possess no transversal degrees of freedom 
we still call them photons. However, thru their interaction with charged fermions 
they may become dynamical fields as exemplified by the Schwinger mechanism. 



Generalized Gauged Thirring Model on Curved Space-Times 3 

is the generally and gauge covariant derivative containing the U{1) gauge 
potential and spin connection, rj^i, = y/—g e^i, denotes the totally antisym- 
metric tensor and TZ the Ricci scalar. The gravitational field g^i, (or rather 
the 2-bein e*^, since the theory contains fermions) is treated as classical 
background field, whereas the 'photons' A^^ 'electrons' V') scalars A and 
pseudoscalars are fully quantized. In two dimensions the electric charge e 
has the dimension of a mass. The other 3 couplings are dimensionless. The 
classical theory is invariant under ^7(1) gauge- and axial transformations 
and correspondingly possesses conserved vector and axial- vector currrents. 

We have chosen this model since it allows to address the above raised 
questions and since it relates to known soluble models for certain values of 
the coupling constants. For example, it contains the gauged Thirring model, 
the Schwinger model in curved space time and the minimal models in con- 
formal field theory as particular limits. Its vacuum structure is non-trivial, 
i.e. it has ^-vacua like more realistic 4-dimensional gauge theories [7]. The 
'photon' aquires a mass = e^/(7r-|- l^r^) via the Schinger mechanism. For 
finite volumes the theory possesses instantons which mimimize the euclidean 
action. These instantons lead to chirality violating vacuum expectation val- 
ues. For example, a non-zero chiral condensate develops which only for high 
temperature and large curvature vanish exponentially. 

For e = all coupling constants are dimensionless and the theory be- 
comes conformally invariant. In this limit the vacuum structure becomes 
trivial. Despite its complexity the general model (1.1) is solvable for arbi- 
trary classical backgrounds ^r^^ and allows for an analytical treatment. This 
in turn enables the entire stress tensor in any curved space, the induced 
currents, their correlators and the equation of state to be constructed. 

The physical role of the coupling constants is the following: The coupling 
of 4> to the transversal current decreases the effective electromagnetic inter- 
action between fermions. For example, the electric charge becomes renormal- 
ized to = e/ ^/Y^\^~g^/27^ , the chiral condensate decreases as ~ (27r-|-5f|)~^ . 
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The mass in the bosonised theory depends on g2. In the ungauged sector 
the Kac- Moody central extension, conformal weights and U{1) charges de- 
pendend on (72. The couphng constant (73 amplifies the Hawking radiation 
which remains thermal for the interacting model. It is (3 + 247Tg^) times as 
strong as that of a free massless scalar field. The central charge and con- 
formal weights depend on ^3. Actually, the weights of the fermionic fields 
become complex for ^ 0. However, g^ does not enter in the finite size 
effects. The coupling constant gi to the longitudinal current weakens the 
long range gauge invariant electron-electron correlators in the one-instanton 
sector (see 6.27). In the ungauged sector it enters in expectation values of 
local operators and in particular in the short distance expansions of the 
fermionic fields and energy momentum tensor. It does not influence the 
thermodynamics of the model. 

Since for particular choices of the coupling constants the model reduces 
to wellknown and wellstudied exactly soluble models there are many earlier 
works which are related to ours. Some of them concentrated more on the 
gauge sector and investigated the renormalization of the electric charge in 
the gauged Thirring model by the four-fermi interaction [8] or the non-trivial 
vacuum structure in the Schwinger model [3,9]. Others concentrated on the 
ungauged conformal sector. Freedman and Pilch calculated the partition 
function of the ungauged Thirring model on arbitrary Riemann surfaces [10] . 
We do not agree with their result and in particular show that there is no 
holomorphic factorization for general fermionic boundary conditions. Also 
we deviate from Destri and deVega [11] which investigated the ungauged 
model on the cylinder with twisted boundary conditions. We shall comment 
on the discrepancies in sections 3 and 7. Other papers which are relevant 
and are dealing with different aspects of certain limiting cases of (1.1) are 
[12], where the thermodynamics of the Thirring model has been studied or 
[5] in which the Hawking radiation has been derived. 
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The paper is organized as follows: In section 2 we analyse the classical 
model to prepare the ground for the quantization. In particular we derive 
the general solution of the field equations, discuss the conservation laws 
and investigate the limiting theories. By employing the graded structure we 
derive the classical Poissson (anti) commutators of the fundamental fields 
with the energy momentum tensor. In the following section we quantize the 
finite temperature model. To avoid infrared problems we assume space to 
be finite. Together with the finite temperature boundary conditions we are 
lead to considering the theory on the 2-dimensional eudidean torus. Due 
to the twists in the fermionic boundary conditions, the non-trivial vacuum 
structure and the associated instantons and fermionic zero-modes the quan- 
tization is rather subtle. Actually we show that some of the results in the 
literature are incorrect. In section 4 the general results are applied to derive 
the partition function of the gauged model. Its dependence on the spatial 
size, temperature and gravitational field is explicitly found. In section 5 we 
show that for equal couplings the gauged model on curved spacetime can be 
bosonized. It turns out that only the non-constant parts of the currents can 
be bosonized and that for this part the wellknown flat spacetime rules need 
just be covariatized. In the following section the chiral symmetry breaking 
is studied. The exact form of the chiral condensate is found. On the flat 
torus the formula simplifles to (6.13). Various limits, e.g. L — > oo, T — > 0, 
T — > oo or ^2 — oo are investigated. By comparing the temperature and 
curvature dependence of the condensate we derive an effective curvature 
induced temperature. In section 7 the thermodynamics of the ungauged 
model is studied. We derive the ground state energy and its dependence 
on the coupling constants, size of the system and boundary conditions. We 
compute the equation of state and our result does not agree with [12]. In 
the last section we investigate the conformal sector of (1.1), that is the 
ungauged model in flat spacetime. Besides the Virasoro algebra the model 
contains an U{1) Kac-Moody algebra. We calculate the important commu- 
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tators and in particular determine the conformal weights and ?7(l)-charges 
of the fundamental fields from first principles. Also we show that the finite 
size effects are in general not proportional to the central charge as has been 
conjectured by Cardy [13]. The appendix A contains our conventions and 
scaling formulae for the various geometrical objects. In appendix B we col- 
lected some useful variational formulae which we have used in this work. In 
appendix C we derive the partition function within the canonical approach. 

2. Classical theory 

2.1 Equations of Motion 

The field equations of the model (1.1) are 

2 V'A = -g^n - g.V^f 

(2.1) 

2 = -g.V^f^ 
V.F^" = er , 

which are the Dirac equation for massless charged fermions propagating in 
a curved space-time and interacting with the scalar and pseudoscalar-fields, 
Klein Gordon type of equations and Maxwell equation. Here j^'^ is the axial 
vector current which is defined by 

f^ = i^ri.i^^vV- (2.2) 
When one decomposes the gauge field as 

A,, = d^a - r]^pdf<f so that Fqi = ^=^VV, (2-3) 

and choses isothermal coordinates for which g^^ = e^^rj^^, then the gener- 
alized Dirac operator reads 

P = eiF-i75G-ia ^ g-iF-i75G+ia^ ^^^^^ 

(2.4) 

F = giX + ea , G = g^cj) -\- e (p . 
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Hence, if ipo{x) solves the free Dirac equation in flat Minkowski space time, 
then 

i;{x) = e'^+'^'^-^''i;o (2.5) 

solves the Dirac equation of the ineracting theory on curved spacetime. The 
vector currents are related as 

= ^^M^ = ^or^oe-^"" = . (2.6) 

The same relation holds for the axial vector current. From -\/— ^V^j'* = 
9nV~~9j^ the conservation of the vector and axial currents follow at once, 

V^j'^ = V^f' = , (2.7) 

expressing the classical ^7(1) x J7a(1) invariance of the model. Using (2.7) 
in (2.1) we conclude that 

2V^A = -g^n and = (2.8) 

or that there is no backreaction from fermions onto scalars. Finally the 
conservation laws imply that the currents are free fields 

V^j^= V^j'^= 0, (2.9) 

which is the reason which accounts for the solubility of the model [14] , even 
in the presence of photons and an external gravitational field. As is well- 
known, for any gauge invariant regularization the axial current possesses an 
anomalous divergence in the quantized model and (2.9) is modified. Thus 
the normal Ua{1) Ward identities in the ungauged Thirring model [8] be- 
come anomalous when the fermions couple to a gauge field. 

In isothermal coordinates the general solution of the field equations can 
be expressed in terms of 6 chiral functions as follows: Introducing lightcone 
coordinates x^ = x^±x^ so that ds^ = e^'^ dx'^dx~ , the solutions of (2.8) read 

X = g-ia + X+{x+) + X-{x~), and (j) = (j)+{x+) + (j)-{x~) (2.10) 
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and depend on 4 chiral functions which are fixed by the initial data on some 
spacehke hypersurface. The solutions of the free Dirac equations depend on 
2 chiral functions as 

*=fM!:|V (2.11) 



In these coordinate system the Maxwell equations (2.1) can easily be inte 
grated and one finds 



a+a_0 = Foi = 2e' 



X X 

j i;liO^_iOd^- j ^i{i)^+{i)d^\ (2.12) 



To go further we must fix the gauge. Conveniently one chooses the Lorentz 
gauge such that q; = in (2.3) and thus in (2.12) determines A^^. We see 
that in isothermal coordinates and this gauge the general solution of (2.1) 
is given by (2.10), (2.12) and (2.5), that is in terms of 6 chiral functions. 

Besides the currents the symmetric energy-momentum tensor of the mat- 
ter fields 

T'^^ = - — (2.13) 

plays an important role in any theory in curved space time. Applying the 
variational identities in Appendix B one obtains after a lengthy but straight- 
forward computation 

+ 2V^0V"0 - ^^"V"(/)Va</> + (</> ^ A) 

- 2g^{gi"'V'^ - V^V'^)A (2-14) 

+ ^i'^(^iV'^A-^2r7"«V,(/)) + iiJi^v) 

where we have introduced the symmetrization A^^^B^^ = \{A^^B^ + A'^B^). 
The first two lines are just the energy momentum of the electromagnetic 
field, charged fermions and free neutral (pseudo-) scalars. The terms con- 
taining second derivatives of A are the improvement terms [15] which are 
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always present when one couples scalars non-minimally to a background cur- 
vature. The remaining terms reflect the interaction between the fermionic 
and auxiliary fields. 

On shell T^^" is conserved as required by general covariance. Using the 
field equations for ij) and A its trace reads 

T% = gln-^F^PF^p. (2.15) 

In particular for (73 = and = it vanishes, and the theory becomes Weyl- 
invariant. As a consequence it reduces to a conformal field theory in the flat 
spacetime limit [16]. It is remarkable that it can be made Weyl invariant 
even when ^37^0. Indeed, without changing the flat spacetime limit we may 
add a nonlocal Wess-Zumino-type term to the action, namely 

S' = S- ^Sp where Sp ^ J ^U^U (2.16) 

the variation of which is 

5Sp = j [Aig^^'n - v^v^^n] + 2^^ (^^7^^ i^h'^ 

The trace of the modifled energy momentum tensor is now zero, and for 
dixu Vnv the Lagrangian corresponds to a conformal fleld theory in 
Minkowski spacetime . 

Choosing the coupling constants appropriately, the model reduces to various 

well known exactly solvable models: 

- For the special choice =5^2 = 6 = and for vanishing gauge fleld the A- 
dependent part of (1.1) is just the Lagrangian of scalar flelds coupled 
to a background charge and and for imaginary g^ describes the minimal 
models of conformal fleld theory [17]. 
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- For ^3 = and g'^ = —g^ = g^ the fermionic sector reduces to the gauged 
version of the Thirring model [18] in curved space time. To see that we 
solve the Klein Gordon equations in (2.1) for the U{1) current which 
yields 

j^ = --d^X--V^''d,cf>. (2.18) 
91 92 ^ 

Inserting this into the Dirac equation we find 

^t'^V^V - yJ'^TmV' = , (2.19) 

which is the field equation of the gauged Thirring model in curved 
spacetime with Lagrangean 

C[i;, ij] = V;z7^V^V - ^rjf. - \f^^F'^'' . (2.20) 

- If we further specialize to = we recover the Schwinger model in 

curved spacetime [9]. 

In the following sections we are lead to consider the euclidean version of 
the model. Then one must replace the lorentzian j'^,gij,v and oj^ by there 
euclidean counterparts. For example, with our conventions (see appendix 
A) the relation (2.2) becomes 

ft^ = -ir^^^j- (2.21) 

and as a consequence the generalized Dirac operator in euclidean spacetime 
becomes 

instead of (2.4). Also, to recover the euclidean Thirring model as particular 
limit of (1.1) we must set gs = and 9^ = gl=g^- 



2.2 Hamiltonicin formalism cind classical conformal structure 
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In this subsection we investigate the Hamiltonian structure of the model 
(1.1) in the conformal limit, i.e. in flat Minkowski space and for vanishing 
gauge field. In the presence of both fermions and bosons it is convenient 
to exploit the graded Poisson structure [19]. We recall, that the equal time 
Poisson brackets are 



){z) ^ 50{z) 5tto{z) 

(2.23a) 

The sum is over all fundamental fields 0{x) in the theory . The sign is 
minus if one or both of the fields A and B are bosonic (even) and it is plus 
if both are fermionic (odd). The momentum densities TToix) conjugate to 
the O-fields are given by functional left-derivatives 

'ts 



ddoO{x) 

A simple calculation yields the following momenta 



= -^V'^ TT^ = g^jl + 2^00 and tta = g^jo + 2^0 A 



(2.236) 



(2.24) 



which form the fundamental Poisson brackets with the fields 

{7Tx{x), \{y)} = 5{x^ - y^), {7r<^(a;), (t){y)] = 5{x^ - y^). 
For the Hamiltonian we obtain 



(2.25) 



iy = y dx^ doijjTr^ + doXTTx + 5o07r<^ - C 



dx^ 77^75^1 V' - igidiX'Ktp^^'ip - ig2di(j)'K^ip + (9iA)^ 
+ {di4>f + ^(tta - ign^i^'ipf + ^(tt^ - ig2T^i,ib'4^f 



(2.26) 



It can be checked that the corrresponding Hamiltonian equations are just 
the field equations (2.1) with flat metric and vanishing gauge potential, as 
required. Since T^^ = (see 2.15) the only non-zero components of T'*^ are 

the lightcone components T_|__|_ and T To continue it is convenient to 

introduce adapted light cone coordinates 
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x^ = x^±x^ so that d± = \{dQ±di), (2.27) 

and the chiral components of the Dirac spinor i)± = i(l±75)V'. Then T__ 
in (2.14) simplifies to 

T__ = - ^(7r^^a_V'+ - a_7r^>+) + 2(a_A)2 + 2(a_0)2 
+ 2^3^^ A + id-{giX + g2(f))'K^+ip+. 

Using the equations of motion one shows explicitly that it is a chiral field, 
i.e. depends only on x~. With (2.28) we can now find the conformal weights 
of the fundamental fields which determine their transformations under in- 
finitesimal conformal symmetry transformations. For that we must calculate 

the commutator of the symmetry generators T/ = J dx~ f{x~)T with the 

fields. The result is 



\ - 

2 



{T/,A} = /a_A-fa_/ 



1 (2.29) 

{Tf, ^+} = fd-^+ + -(1 - ig^gs)^+d-f 

{Tf, = fd-iPl + ^(1 + ig,gs)iPld_f. 

Whereas and ■0+ are primary fields, A is not. Actually, the non-primary 
character of A is very much linked with the gfs-dependent term in the 
transformation of the Dirac field. To see that more clearly we note that 
under an infinitesimal left conformal transformation generated by Tf = 
f dx~^ f{x~^)T++ the scalar and fermi field transform as 

{ff,X} = fd_X-^d-f and {ff,,p+} = fd-,p+-ig^gsiP+d-f. (2.30) 

Since ■0+ is not any longer a scalar under left transformation the term 

dx'^dx~ {2itl}\{dj^ — igid+X)ip+^ 
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appearing in the action is only conformally invariant because A transforms 
inhomogenously like a spin connection. It may be surprising that the sym- 
metry transformations depend on the coupling constant which is not 
present in the flat space time Lagrangean. However, the same happens for 
example in 4 dimensions if one couples a scalar field conformally, that is 
non-minimally, to gravity. Although the Lagrangeans for the minimally and 
conformally coupled particles are the same on Minkowski spacetime, their 
energy momentum tensors are not. The same happens for the conformally 
invariant nonabelian Toda theories wich admit several energy momentum 
tensors and hence several conformal structures [20] . 
The current transforms as 

{T/,j_} = /c»_j_+j_c>_/ (2.31) 

and the energy momentum tensor as 

{T^,T__} = fd.T__ + 2T__5_/ - gldlf. (2.32) 

Recalling that a primary field O with weight h transforms as 

{Tf,0} = fd_0 + hOd_f 

and comparing with the above results we have found the following structure: 

- The pseudoscalar field is primary with /j,,^ = 0. The scalar field A is 
only primary for ^3 = in which case h\ = Q. 

- The Dirac field ijj^ is primary with h^j^ ~\{^ ~ WiQs)- The conformal 
weight is real for imaginary g^. 

- The current is primary with weight 1. 

- Already on the classical level the energy momentum tensor is only quasi- 
primary. The corresponding Virasoro algebra (2.32) has central charge 
c = 247r^|. 
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3. Quantization of the generalized gauged Thirring 



In this section we quantize the general model (1.1) in curved space-times. 
The results are then applied in the following sections, where we calculate 
the partition function, ground state energy, equation of state and certain 
correlators of interest and their dependence on the chemical potential, vol- 
ume of space, temperature and background metric. To do that we couple 
the conserved U(l)-charge to a chemical potential We enclose the system 
in a box with length L to avoid infrared divergences. To investigate the 
temperature dependence the time is taken to be purely imaginary in the 
functional approach [21]. The imaginary time varies then from zero to 
the inverse temperature /? and we must impose periodic- and antiperiodic 
boundary conditions for the bosonic- and fermionic fields, respectively. Thus 
to study the finite temperature model we must assume that space-time is 
an euclidean torus [0,(3] x [0,L]. 

To see how the partition function and correlators depend on the gravi- 
tational field we assume that the torus is equipped with an arbitrary met- 
ric with euclidean signature or equivalently with a 2-bein e^a- The curved 
gamma matrices are 7^ = e^al"' and in particular 75 = —^r]^^^'^^'^^ — as 
is constant (see appendix A for our conventions). We can always choose 
(quasi) isothermal coordinates and a Lorentz frame such that 



where t = ti + iro is the Teichmueller parameter and a the gravitational 
Liouville field. Space-time is then a square of length L and has volume 



model 




or 




(3.1) 



V = Jq dPx^Jg. We allow for the general twisted boundary conditions for 



the fermions 
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The parameters ctj and /?j represent vectorial and chiral twists, respectively. 
We could allow for twisted boundary conditions for the (pseudo) scalars as 
well, e.g. (p{x^ + nL,x^+mL) = cl){x^,x^) + 2n{m+n). However, to recover 
the Thirring model for certain values of the couplings we assume that these 
fields are periodic. For a = 0, t = iP/L and ao = Pq = the partition 
function has then the usual thermodynamical interpretation. Its logarithm 
is proportional to the free energy at temperature T=l/(3. 

3.1 Fermionic path integral 

Twisted boundary conditions as in (3.2) require some care in the fermionic 
path integral. Indeed the fermionic determinant is not uniquely defined when 
one allows for such twists. The ambiguities are not related to the unavoidable 
ultra- violett divergences but to the transition from Minkowski- to Euclidean 
space-time. To see that more clearly let denote the set of fermionic fields 
in Minkowski space-time with chirality ±1. Since both the commutation 
relations and the action do not connect and S~ we can consistently 
impose different boundary conditions on and S~ . On the other hand, in 
the euclidean pathintegral for the generating functional 



exchanges the two chiral components of '0, i.e. p : . Thus, in con- 

trast to the situation in Minkowski space the two chiral sectors are related 
in the action. Of course, the eigenvalue problem for ip is then not well 
defined. This is the origin of the ambiguity in the definition of the deter- 
minant. It is related to the ambiguities one encounters when one quantizes 
chiral fermions [22]. To solve this problem we shall analytically continue 
the well-defined determinants in the untwisted sector f3 = to /? 7^ 0. The 




(3.3) 



the Dirac operator 
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resulting determinants do not factorize into (anti-) holomorphic pieces and 
differ from previous ones in the literature [10]. In appendix C we give fur- 
ther arguments in favour of our result by calculating the determinants in a 
different way. 

Let us now study the generating functional for fermions in an external 
gravitational and gauge field and coupled to the auxiliary fields. For that 
we observe that on the torus the decompositon (2.2) of the gauge potential 
generalises to 

An = K+ X^'" ^ ^""^ ~ Vn^d'^'fi, (3.4a) 

where the last 3 terms are recognized as Hodge decomposition of the single 
valued part of A in a given topological sector, that is the harmonic-, exact- 
and coexact pieces. In arbitrary coordinates the toron field obeys the 
harmonic conditions V^t'* = t[^^i/] =0. It follows then that in isothermal co- 
ordinates must be constant. The role of the toron fields has recently been 
emphasized within the canonical approach [23]. In the Hamiltonian formu- 
lation they are quantum mechanical degrees of freedom wich are needed for 
an understanding of the infrared sector in gauge theories. Also, in [24] it 
has been demonstrated that the Zjv-phases of hot pure Yang-Mills theories 
[25] should correspond to the same physical state if one takes care of the 
toron fields. 

The first term in (3.4a) is an instanton potential which gives rise to a 
nonvanishing quantized fiux or integer-valued instanton number k: 



= e J Foi = e J E = e J = 27rk. (3.5) 



As representative in the A;-instanton sector we choose the, up to gauge trans- 
formations, unique absolute minimum of the Maxwell action in (1.1). It has 
field strength eE^ = ^Jg^jV . As instanton potential we choose 

eA^ = ei^ -$r]/d^X, where ei^ = -^{x\o) (3.46) 
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is the instanton potential on the flat torus with the same flux but field 
strength yfg^jV . The function x is then determined (up to a constant) by 

- V^J = V?Ax. (3.4c) 

The solution of this equation is given by 

X{x) = --i(le-2-)(x) = iy d2yV^Go(x,y)e-2-(^), (3.4d) 

where 

G„(.,»)=(.|-L|,)=;^«£MM (3.6a) 

A„>0 " 

is the Greenfunction belonging to —A. In deriving (3.4d) wc have used that 
;g-(<?/y) =0 which follows from the spectral resolution (3.6a) for the Green 
function in which the constant zero mode 0o = 1/ of A is missing. 
Note that 2-dimensional gauge theories are not scale or Weyl invariant as 
4-dimensional ones are. For that reason the instantons on conformally flat 
spacetimes are not just the 'flat' instantons. 

To be more explicit we relate Go to the Greenfunction Go on the flat 
torus with the hatted metric [26] 

Go(a;,y) = -— log|— -[j ^1(0,T)|^ where i = x-y. (3.66) 

For that we note that due to the missing zero-mode in (3.6a) the usual flat 
spacetime equations for the Greenfunctions are modifled to 

-A^Go(a;,y) = — , -A^Go(a;,y) = y= (3.7a) 

Furthermore one sees at once that both Green functions annihilate the cor- 
responding constant zeromodes 

j (fy.My)Go{x,y) = J (fy^gGo{x,y) = 0. (3.76) 

From these two equations one concludes that Greenfunction on the curved 
torus is related to the flat one (3.6b) as 
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Go{x,y) ^ Goix.y) + j (fu<fvsJg{u)g{v)GQ{u,v) 

-yj Go{x,u)-\/g{u)d'^u - y j (fu-\/g{u)Go{u,y) 



(3.8) 



and this replaces the infinite space relations Gq = Gq [27]. 

Our choice for the instanton potential (3.4) corresponds to a particular 
trivialization of the t/"(l)-bundle over the torus [3]. In other words, the 
gauge potentials and fermion fields at {x^, x^) and {x^, x^+L) are necessarily- 
related by a nontrivial gauge transformation with windings 
^^(^o, x^ + L)- yl^(a;0, x^) = d^a{x) 

{3.9a) 

For the choice (3.4b) we find 

ea{x) = ~x^. (3.96) 
Ij 

Note that A is still periodic in x^ with period L and ij) still obeys the first 
boundary condition in (3.2). Our trivialization differs from the one chosen 
in [28] and so do our instantons and fermionic zero modes. 

Similarly as for the gauge potential we must add a harmonic piece to 
the aiixiliary vector field to which the fermions couple in (1.1), so that 

27r 

= —goh^ + ^i^^A - g2ritiud''(t> (3.10) 

appears in the Dirac operator in (1.1) on the torus. A and couple to the 
divergence of the vector and axial vector currents. The harmonic fields 
couple to the harmonic part of the current and are needed to recover the 
Thirring model in the limit gQ=gi=g2- Also, we shall see that and /i^ are 
essential to obtain the correct answer for the thermodynamic potential. Note 
that contains no instanton part since it couples to the gauge invariant 
fermionic current. 

Finally we introduce a chemical potential for the conserved U{1) charge. 
In the euclidean functional approach this is equivalent to coupling the 
fermions to a constant imaginary gauge potential Aq [29] . 



Generalized Gauged Thirring Model on Curved Space-Times 19 

Inserting the above decompositions and the chemical potential into the 
Dirac operator finally yields in isothermal coordinates 

p = = e'^+i^^G+^^'^-i^ p e-^F+^,{G+'Px)+h'^ ^ ^here 

P = l^{d^. + - ieA^ -—[H^ + //^]), (3.11) 
H^ = et^+ goh^ and //^ = -^-^A* ^^o- 

Here u> is the spin connection belonging to e^a- It vanishes for our choice of 
the reference zweibein. is the instanton potential (3.4b) on the flat torus. 
The scalar and pseudoscalar functions F, G and x have been introduced 
in (2.4) and (3.4d). In the chosen coordinates t and h and hence H are all 
constant. In [3] it has been shown that p possesses zero-modes of definite 
chirality and their chirality is given by the sign of k. They are crucial in any 
correct quantization. For example, if one would leave out instanton sectors 
in which ip has zero-modes then the cluster property would be violated. 

In a first step wc quantize the fermions in the fiat instanton and har- 
monic background and reference metric g^^, that is we assume p ^ p in 
(3.3). The dependence on the remaining fields F,G,x ^^d cr, that is the 
relation between Zp and Zp, is then found by integrating the chiral and 
trace anomalies [30] and exploiting the relation (3.11) between p and p. 
We expand the fermionic field in a orthonormal basis of the Hilbert space 

^(^) = ^an1pn+{x) + ^bn1pn-{x) 

n n 

V'Ua;) = ^ anxi+ i^) + Kxl- (x) , 

n n 

where an,bn,anjbn ^tre independent Grassmann variables. 



Topologically trivial sector 

For k — or vanishing instanton potential we can immediately write down 
a basis 

^^^{x)^-^e'^P-'^^e±, where (p±), = — (- + a, ± A + n^), (3.13) 



(3.12) 
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and e± are the eigenvectors of 75. The i(^n+ ^ind V'n- must obey the S'^ 
and S~ boundary conditions, respectively. These boundary conditions fix 

the admissable momenta in (3.13). Since the Dirac operator maps 
into the Xn± must then obey the same boundary conditions as the 'i/'n=F- 
Thus 

Xn±{x) is obtained from ipn±{x) by exchanging and p^ . It foUows 
then that 

ipipni = A^Xn=F (3.14a) 



with 



K =^[^{l + ai+(3i+ m) -{lr + ao + (3o + no)] 
tqL z z 



(3.146) 



Here we have introduced = a^— i^^— Substituting (3.12,14) into the 
generating functional (3.3) and applying the standart Grassmann integra- 
tion rules we arrive at 

ZF[v,fj] = det ip e-/*^(")^("'^)^(^\ detip ^Hx+X-, 

n 

A. X \^ ^ i^n+ix)xi-iy) , V^n-(a^)Xn+(^) ^ 

S{x, y) = 2^ ( -+ + ). 

S is the fermionic Green function in the 0-instanton sector. Note that both 
the 'eigenvalues' and the Green function depend on the Teichmueller pa- 
rameter, harmonic potentials, twists and chemical potential. 

We proceed to calculate the infinte product or generalized determinant 
in (3.15). This is one of the central points of our article and for non-zero 
chiral twists and chemical potential our result deviates from previous ones 
[10]. Actually the twists and chemical potential are related as one can see 
from (3.14). One may be tempted so identify 



det(i:'+i:'_) - Yl X+X~ and det D+ det i:>_ ~ A+ J]^ A;;, (3.16) 
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and thus conclude that the determinant is a product, /(t)/(t), that is 
factorizes into holomorphic and anti-holomorphic pieces (the overall factor ~ 
I/tqL in the eigenvalues (3.14b) drops in the infinite product, since the torus 
has vanishing Euler number). However, the infinite product in (3.15) must 
be regularized and the two expressions in (3.16) may difli^er. In conformal 
field theory [26] one is naturally lead to consider the individual chiral sectors 
and thus finds holomorphic factorization. For Dirac fermions one uses p'^ 
to regularize the product and this leads to the determinant of the product 

To continue we recast the infinite product in the form 

nA+A;= n (^)V^(^ + C;. + n^)(^ + c. + n.) (3.17a) 

where g^'^ is the inverse of the reference metric (3.1) and 

c^ = a^ + ifj^''p,, where (fy^ '^) = -1 (^^^^ -rf ) ' ^^'^^^^ 
The point is that for real c^, that is for vanishing chiral twists and 
chemical potential (see the definitions of below (3.14b) and /x^ in (3.11)) 
the zeta function defined by 

C(«) = E(^nA;)"^ (3.17c) 

n 

has a well defined analytic continuation to s < 1 via a Poisson resummation. 
An explicit calculation yields [3,31,38] 

detm ^ {l[>^+X-),,g = e-^'(^)l^=°, where 

(3.18) 

L|?7(t)|^ L Co -I L Co J J 

However, for complex c^ the Poisson resummation is not applicable and 

C'(0) cannot be calculated by direct means. To circumvent these difficulties 

we note that the infinite product (3.17c) defining the C- function for s>l is 

a meromorphic function in c. Thus we may first continue to s < 1 for real c^ 

and then continue the result to complex values. Using the transformation 

properties of theta functions the resulting determinant can be written as 
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det(z^) 



:0 



-ai + Pi 
ao - Po 



(O,t)0 



-ai - /3i 
ao + Po 



(0,r). 



(3.19) 



It can be shown that this determinant is gauge invariant, i.e. invariant un- 
der Q!^+l, but not invariant under chiral transformatins, — > P^+l, 
as expected. Furthermore it transforms covariantly under modular transfor- 
mations T — > r -|- 1 and r — I/t- In other words, det ip is invariant under 
modular transformations if at the same time the boundary conditions are 
transformed accordingly. The exponential prefactor is needed for modular 
covariance and is not present in the literature [10]. It correlates the two 
chiral sectors and will have important consequences. In the appendix C we 
confirm (3.19) with operator methods. 



Topologically nontrivial sectors 



For definiteness we assume k > 0. Then ip possesses k zero-modes V'o+j P — 
1, . . . ,k with positive chirality and boundary conditions. Together with 
the excited modes 'ipn+ they form a basis of S'^ . Thus we must add the zero- 
mode contribution ^ CnipQ.^. to ip in (3.12). Similarly we must add c„Xo+ to 
ifj^ in (3.12). The zero modes Xo+ i^i V'^ must obey S- boundary conditions 
(see the discussion below (3.13)). Thus the zero- modes in the expansions of 
ijj and t/j^ have the same chirality but obey different boundary conditions. 
This is required for the zero- and excited modes to form a complete basis 
and is consistent since ip does not relate the zero-mode sector of S'^ with 
S~. The Grassmann integral over the variables belonging to the excited 
modes is performed as in the trivial sector. Also, the integration over the 
Cn and Cn can easily be done and one obtains 



Zf[v,v] 
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k 

"''o+AAo-l-' 



det'ip= n A+A-, Mx,y)= ^ ftll±^^l^ + i+^-)). 

(3.20) 

Note that the excited Green function anticommutes with 75. 
To calculate the determinant we observe that 

= ( "f' dX ) = ^^^-V^S""^^ - \^ + l^-^'.-T. (3.21) 
simplifies in the instanton background and on the flat torus to 

= -g^'^D^D, - t^,. (3.22) 

In other words, it is the same in the two chiral sectors, up to the constant 
This observation allows one to reconstruct the spectrum of —p"^ as 

follows: 

First note that we can define two sets of normalizable zero-modes in the posi- 
tive chirality (75 = 1) sector. One containing k modes obeying the bound- 
ary conditions and the other consisting of k modes with the S~ boundary 
conditions. The first set is admissable and are just the k zero modes -00-1- 
appearing in (3.20). The other k zeromodes Xq_^_ appear also in (3.20) but 
the Dirac operator does not act on them. But because of (3.22) they are 
at the same time eigenmodes of —p'^ in the negative chirality sector with 
the correct S~ boundary conditions and eigenvalues 2$/V. Of course, the 
t/'q^. are also eigenmodes in the 75 = —1 sector but with the wrong boundary 
conditions. However, applying l)_ to them produces k eigenmodes in the 
positive chirality sector with the correct boundary conditions and eigenval- 
ues 2<P/V. This procedure may now be iterated and one ends up with the 
following spectrum of —p'^: 

^2 ^ f ^ deceneracy = k 

1 2n^/V deceneracy = 2k. \ ■ > 



24 A. Dettki, I. Sachs and A. Wipf 

With the exphcit spectrum at hand we can compute the zero-mode trun- 
cated determinant with zeta-function methods and find [3] 

dei'm = (— ) . (3.24) 

We proceed with computing the zero modes of p^. For that we note that 
the operator commutes with the time translations which leads to the ansatz 

~^^^2.ic,x°/L^2.iH,.^/L^^^^i^^^^ Cp=p+(^ + ao + /3o), (3.25a) 

where we have assumed A; > 0. The choice of Cp is dictated by the time- 
like boundary conditions in (3.2). Inserting this ansatz into the zero mode 
equation P'^Xp — yields 

where y = x'^ + ^{cp - Hq - jio). 

This is just the differential equation for the ground state of a generalized 
harmonic oscillator to which it reduces for r = ztq. The solution is given by 

r ^ r 1 ^/ TT N-i2 
ip = exp ^ - ^T^[x +-{Cp-Ho- /io) I 

These functions do not obey the boundary condition (3.9), but the cor- 
rect eigenmodes can be constructed as superpositions of them. For that we 

observe that 

so that the sums 

(2fcTo)^ ^;g o^i^Ho-ao-h)f3i 



= / e 



(3.256) 



^ -2i7r(n+p/fc)(|+ai+/9i-Hi), 

nez 

where p = 1, . . . , /c, obey the boundary conditions and thus are the k required 
zero-modes. We have chosen the phase such that the accompanying zero- 
modes in (3.20) are just 
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Xl+{x,p,...)=i,^4x,-p,...). (3.25c) 

Actually the product Xo+V'o-i- only determined up to a (possibly a and 
/^-dependent) phase. But gauge invariance requires that 

X 

Xo+{x) il^o+{x) and Xo+(^) [^^ J ^] '^o+iv) 

y 

are both invariant under — > a^ + n^ and et^ — > et^^ + n^, where the 
are integers. This almost fixes the relative phases of the zero-modes in 
(3.25b) and (3.25c). Also, the overall factor normalizes these functions to 
one. Modes with different p are orthogonal to each other, so that the system 
(3.25b) forms an orthonormal basis of the zero-mode subspace. For /c < 
the zero-modes are the same if one replaces jS^ by —(3^ and e+ by e_. 

Integrating the chiral and trace anomalies 

To relate the determinants of ip and ilp we introduce the one-parameter 
family of Diracoperators 

^ gT[iF+75(G+^x)-i(T] p ^T[-iF+'i^{G+^x)+h<^] (-3 26) 

which intepolates between p and p [30] . The r-derivative of the correspond- 
ing determiants is determined by the chiral and trace anomaly. An explicit 
calculation yields 

1 

l°g^^ = 7- /^^ / V^tral(27.[G + ^x]-^)+logdet:^. (3.27) 

^ 

Here g'^ is the determinant of the deformed metric gl^^ = e^'^'^g^i,, and 



1 . . „ 1 



12 



a[ = --i?^ + 75TA-G+^ (l-T)V^^+rV^- 75 ( 3.28) 



is the relevant Seeley-deWitt coefiicient of p'^. Furthermore, is the norm- 
matrix of the zero-modes V'o+ (3.25b). Since those are orhonormal it is 
just the /c-dimensional identity matrix. J\f^ is the norm-matrix of the zero- 
modes of ip which are related to the V^g-i- ^ 
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ijjP^ = eiF-75(G+^x)-^a^p_^ (329) 

as follows from (3.11). Inserting (3.28) into (3.27) one finds the following 
formula for the determinant in arbitrary background gravitational and gauge 
fields: 

det'i^ = det ^ det'(i^) exp (^^) 

"... (3-30a) 



•exp / V^GAG+ I j^9G^^j V^x), 



where 



Sl = j ^/^[na-aka] (3.306) 

is the Liouville action. In deriving this result we used that J ^/gx — ^- Ac- 
tually, for our reference metric the Ricci scalar TZ vanishes and the Liouville 
action simplifies to — J y/gaAa. However, as it stands the formula (3.30) 
holds for arbitrary reference metrics and arbitrary Riemannian surfaces. 

As expected for a gauge-invariant regularisation, the function F and 
thus the pure gauge part of the vector potential does not appear in the 
determinant. 

For later use we also give the analogous formula for the zero-mode trun- 
cated scalar determinant [32] 

■ Va ( 1 

'V' 

This completes the computations of the determinants. 

The generating functional for the full theory is then obtained as follows: 
First one notes that the formulae (3.15) and (3.20) for the fermionic func- 
tionals still hold without hats. Thus to calculate the functionals in arbitrary 
gauge-, auxiliary- and gauge fields we need to know the Greenfunctions, de- 
terminants and zero-modes in these backgrounds. 

To relate the fermionic Greenfunctions S in the different topological sectors 
to the hatted ones we define 

Si{x,y)^e-^^^^ S{x,y)e-~^^y\ g = -iF + ^,{G + ^x) + (3.32a) 



det'^ (- A) = det'^ (- A) (|) ^ exp ( - ^S^) . (3.31) 
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On the infinite space we would have S = Si [27] . However, if the Diracop- 
erator possesses zero modes this simple relation is modified to to 



S{x,y) = Si{x,y) + / Po{x,u)Sl{u,v)Po{v,y)^y g{u)g{v)(fud 



2„ 



V 



(3.326) 

Si{x,u)Po{u,y)y/g{u)d'^u - I Po{x, u)Si{u, y)^g{u)(fu, 



and this formula should be compared with the analogues one for scalars 
(3.8). Here Pq is the orthonormal projector onto the zero modes. For gauge 
fields with vanishing flux S — Si. Together with the relation (3.30) between 
the full and hatted determinant and the explicit form (3.18,19) for deti^ 
this yields the fermionic generating functional in the various toplogical sec- 
tors. 

In the trivial sector one finds explicitly 



1 



T 



-Cl 



ZFhfi] = j::j^e / (o,t)0 ,Mo,t) 



-Cl 



J ri{x)S{x,y)ri{y) _ 



/I 1 r ^ 



By using the scaling properties of the Ricci-scalar and Laplacian (see ap- 
pendix B) the exponent can be rewritten as 



I67r J A 27T J 



967r 

which makes clear that the resulting functional is diffeomorphism invariant. 
Here we used that TZ integrates to zero or that the Euler number of the 
torus vanishes. On the sphere or higher genus surfaces the last formula is 
modified. 

To relate the hatted and full functionals in the non-trivial sectors one recalls 
that in the formula (3.20) for the full partition function (without hats) 
one must use orthonormal zero-modes. These can be expanded in terms of 
the un- normalized modes V'o-i- Xo-i- defined in (3.29). Inserting these 
expansions into (3.20) yields the inverse square roots of the determinants of 
the corresponding normmatrices jV^ and A/"^ which partly chancel det J\f^ 
in (3.20). Thus one ends up with 
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where the and Xo+ ai'e the un-normahzed zero-modes (3.25). 
3.2 Bosonic path integral 

To arrive at the generating functional for the complete theory we must 
finally quantize the photon and auxiliary fields and (see (3.10)). 
For that we insert the decomposition (3.4a) into the bosonic part of the 
(euclidean) action (1.1). This results in 

r n X ^3.35) 

+ / y^(^-v.AV-AAA-0A(/.-^3i?Aj. 

The term quadratic in the h field is not present in the action (1.1) on 
Minkowski space-time. But on the torus h is part of the Hodge decompo- 
sition of and thus on the same footing as dX and rjdcj). Since Sb and 
the fermionic determinants are both gauge invariant and thus independent 
of the pure gauge mode a in (3.4a), it is natural to change variables from 
to {(p,a,tn,<P) in each topological sector. One can show [3] that this 
transformation is one to one, provided 

JVgV=JV9(^ = and et^e[0,l]. (3.36) 

The measures are related as 

VAf, = dtodtiV<fVa, where J = (27r)2det'(-A). (3.37) 

k 

The Jacobian J is independent of the dynamical fields. In expectation values 
of gauge invariant and thus -independent operators the a-integration cancels 
against the normalization. This is of course related to the fact that in QED 
the ghosts decouple in the Lorentz gauge. 
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Finally observe that via the derivative couplings to the fermionic current 
[33] we introduced artificial degrees of freedom. The relation between in 
(3.10) and the fields (0, A, h^) is only one to one if we impose the conditions 
similar to (3.36), namely 

4'=^ J \/94^ 0: A=0 and /i^ G [— oo, oo]. (3.38) 

There is no restriction on the harmonic part of the auxiliary field, since 
is not a gauge field. The constraints are imposed in the functional integral 
as 

J dhodhiV(f>VX6{(t))6{X) (3.39) 

The normalization by the volume in (3.38) is needed such that the con- 
straints and hence the partition function are both dimensionless. For exam- 
ple, expanding in eigenmodes of the Laplacian as 

= ao0o + y2 ^ri(pn , where (po = 

n>0 

is the zeromode, one finds the dimensionless partition function 

/ V4> 5{4>) et'^'t' = W ^} (3.40) 
J det'2(-A) 

for free bosons. 

Constraining the mean field to zero as in (3.40) is equivalent to fixing the 
field at an arbitrary point ^ on the torus to zero [34] 

j Vcf>5i4>)--- = j V4 5{m)---- (3.41) 
This can be seen as follows: 



/ 



V(j> 5{m) ■■■ = j du5{4>-u)V<j> 5{m) ■■■■ (3.42a) 



Now one shifts the field as Using that the action is left invariant 

by this shift, the measure becomes 

j duVcP 5{(p)5{(p{C) + u) . . . = Jv(P ■ ■ ■ (3.426) 
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which shows that the two constraints are the same. When integrating over 
the airxihary fields it is always understood that the divergent zeromodes are 
suppressed as in (3.39). 

4. Partition function of the generalized Thirring 
model 



As a first application of our general results we calculate the partition func- 
tion of the theory (1.1). To compute it we must put the sources rj and ry in 
(3.3) to zero. Then it is evident from (3.34) that the non-trivial sectors do 
not contribute and hence we may assume <? = 0. Thus the partition function 
is given by 

Zo = J J (ft(fhV(pV(l)VX Zf[0,0] e-^«[*=o], (4.1) 

where J is the Jacobian of the transformation (3.37). Zp is the fermionic 
partition function (3.33) in the trivial sector and the integration is over fields 
obeying the conditions (3.36,38). Now we perform the various integrals in 
turn. 

Integration over the harmonics: 

By using the series representation of the theta functions one computes 



1 






Since the result appears always together with the ry-function factor in (3.33) 
it is convenient to introduce 

_ 1 1 
V2^|r/(r)P 

in the following expressions. The result (4.2) does not depend on the /i-field 
and hence the /i- integration in (4.1) becomes Gaussian. It yields a factor 
l/47r so that 
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Zo = 7r«det'(-A) J Vs{ip<l>X) e^/v^«^G-5s[h=o]^ ^^ ^^ 

where G has been defined in (2.4). We inserted the exphcit expression (3.37) 
for the Jacobian. It is interesting to note that aheady the toron-integration 
in (4.2) washes out the dependence on the boundary conditions and chemical 
potential. We shall comment on this point later on. 
Integration over X and 0; 

The integral over A, subject to the condition (3.38), modifies the Liouville 
factor and yields one inverse square-root of the determinant of — 2A in (4.3). 
To continue we recall the scaling formula for the determinant of A [35] : 

det'f— aA) , ^, , , r 1 /" i /, , x 

log = log a • C(0) = log a • [— J ai-p], (4.4a) 

where p is the number of zero modes of the operator. On the torus J ai = 
and we find 

det'( - aA) = - det'(-A). (4.46) 

a 

Using this scaling property the X-integral together with (3.40) we obtain 
Zo = ^«7rdet'^(-A) e'^9i+i/24-)s, J J ^GAG-Sslh^x^o] ^ 

(4.5) 

To quantize the field we need to recall that G = e(p + g2(f). Since 
(pAif ~ (A^jA^), the anomalous term ~ J GAG in the exponent con- 
tains an explicit photon mass term with bare-mass e/y^. However, when 
quantizing the (p field this mass is renormalized. This can be seen explicitly 
in the resulting expression for the partition function after the (^-integration 
has been performed 

J 

where the renormalized photon mass is 

-^2 27r 



ml = — ^ ^ . (4.7) 

^ TT 27r + ^2 ^ ^ 



Integration over ip: 
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The zeta-function regulated determinant which one obtains when perform- 
ing the integral (4.6) factorizes 

det'(A2 - m^A) = det'(-A) • det'(-A + m^). 

This factorization property is not obvious since all determinants must be 
regulated. But it holds for commuting operators and in the zeta-function 
scheme. Then the partition function simplifies to 

= ^^^^(det'(-A)det'(-A + m^))-^ exp ((^3' + 3^)^^) • (4.8) 

We can go further by using (3.31) and the known result for the determinant 
of A [26] which together yield 

det'^(-A) = roL|ry(r)pyi exp ( - ^S^) (4.9) 

which finally leads to 

Zo = V2^— r^--rj 1 Gxp({-^ +gi)SL] (4.10) 

m^ro|r7(T)|4 det'^(-A + m2) "^^^I^tt ) ^ > 

for the partition function of the general model (1.1) on curved spaces. It 
shows explicitly that in the topologically trivial sector the theory should be 
equivalent to a theory of free massless and massive bosons with mass m^. 
It is interesting to follow the various contributions to the explicit dependence 
on the gravitational field since they contribute to the Hawking radiation. For 
that we recall that when one quantizes a conformal field theory with central 
charge c in an external gravitational fields one ends up with the Liouville 
term, Z ~ exp[c5'L/247r] [32]. Thus the fermions contribute with c—1, as 
expected. The (j) and A field contribute with 1 and 1 -|- 2A7rg^, respectively. 
However, the Jacobian combined with the conformal part of the gauge sector 
contribute with c = — 1 and we are left with a total central charge c = 
2 + 247rfir|. Of course, the gauged model is not conformally invariant and 
the breaking is manifest in the massive determinant in (4.10). The partition 
function of the ungauged theory is (4.6) multiplied by an inverse determinant 



Generalized Gauged Thirring Model on Curved Space-Times 33 

(the missing Jacobian) and without (^-integration. In this hmit one obtains 
a conformal theory with central charge c=3 + IAttq^,. 

By using an elegant result of Christensen and Fulling [36], that relates 
the conformal anomaly to the asymptotic Hawking flux, one concludes that 
the Hawking radiation of the ungauged model is 3 + 2477^1 times that of 
free massless scalars. For the gauged model the Hawking radiation is still 
thermal and consists of massless and massive particles. 

The appearance of in (4.6) should be interpreted as renormalization of 
the electric charge induced by the interaction of the auxiliary fields with 
the fermions. After summing over all fermion-loops this leads to an effective 
coupling between the photons and the (/)-field and in turn to a modified 
effective mass for the photons in (4.6). In the limit g2 ^ ^ this mass tends 
to the well-known Schwinger model result, — > ej sfii [37]. 

We have already mentioned that the chemical potential coupled to the elec- 
tric charge has completely disappeard from the partition function. This does 
not come as a surprise since the only particle in the gauged Thirring model 
is a massive neutral boson. This has no charge which may couple to the 
chemical potential. If the partition function would depend on then the 
expectation value of the charge would not vanish, in contradiction to the 
integrated Gauss law. 

Also note that the final result is independent of the chiral and non-chiral 
twists. The normal twists have been wiped out by the toron integration. 
Since the chiral twists are equivalent to chemical potentials the partition 
function should not depend on them in order not to violate gauge invariance. 
If one would assume holomorphic factorization for the fermionic determinant 
[11] then the partition function would depend on the chiral twists. 

We conclude this subsection with deriving an explicit formula for the parti- 
tion function on the fiat torus. Applying the results in [38] one obtains for 
the massive determinant 
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det'(-A + m^)^ = — e-^^'(°), (4.11a) 
^ rriy 

with 

C'(0) = E :^^^,Mm,L./^)) - (4.116) 

where (n, n) = gijffn^ is the inner product taken with the reference metric, 
and the sum is over aU (n*) e Z2 with the origin excluded. For g^i, = 
Snui in which case the partition function has the usual thermodynamical 
interpretation, the result reduces to one derived previously by Ambjorn 
[39]. In addition, if L approaches infinity we recover a result in [29]. The 
free energy for ti = and on flat space simplifles then to 

i^=-^logZ=^C'(0). (4.12) 
with C'(0) from (4.11b) and the particular choice for the parameters. 

5. Bosonization of the gauged Thirring model 

In the classical analysis we have already seen that in the limiting case ^3 = 
and gi = g2=g the general model reduces to the gauged Thirring model. In 
this section we show that the same is true for the quantized theory on the 
torus if in addition we set go = g- More precisely, the Hubbard-Stratonovich 
transform [40] of the Thirring model is just the derivative coupling model 
(1.1) with identical couplings. In the process of showing that we shall arrive 
at the bosonization formulae for the gauged Thirring model on the curved 
torus. We shall see that only the non- harmonic part of the fermion current 
can naively be bosonized and that for this part the rules of the ungauged 
model on flat space time [41] need just be covariantized. 

For that we calculate the partition function (4.1a) in a different order. 
First we integrate out the auxiliary flelds. In order to understand the role 
of A and we introduce sources for them. Thus we study the generating 
functional for the correlators of the auxiliary flelds 
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Z[^X]= J V{X(^hijA^)e-^+I M^^+^'t'] . (5.1a) 

Here 

S=-i j \f9i'''P^ + Sb{9z^^\ (5.16) 

is the action of the fuU theory. Tj) is the Dirac operator in (3.11) with aU 
couphngs set equal and Sb the bosonic action (3.35). Since A and integrate 
to zero (see (3.38)) we may assume the same property to hold for the sources. 
Also, since there are no fermionic sources only configurations in the trivial 
sector contribute, so that there is not instanton potential in (3.11) and 
hence <? = in (3.35). The integration over the auxiliary fields is gaussian 
and yields 

Z=M) y V{i,A^)e-^- exp ^ V^[-i(^le+C^C) + |(^^jrM+C^i5^) 

(5.2) 



where 

St — — 7 



- J ^(f^.F^'^ - li^^pi, - ^rj^^ (5.3a) 



4 

is the action of the gauged Thirring model on curved space-time and 
comes from the integration over the auxiliary fields. 

Let us first consider the partition function, that is set the sources to zero. 
Comparing (5.2) with (4.6) and using (4.9) we easily find 



P(V't)e-^'^ = \ll + ^^~^-^ ^""^^"^^ j ^7 '^(7) , (5.4a) 

where 7 is the mean field (see 3.38) and we used (3.37) and (4.4b). The 
action for the neutral scalar field 7 is found to be 

S,= i/v5^'..a'-T-^^7jJp7^/v57A.. (5.46) 

Since (5.4) holds for any </? (and thus for the non-harmonic part of any A^, 
because of gauge-invariance) we read off the following bosonization rules: 
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1 



^I'^d^-f and j'i 



1 



3 



0r yrT^2/^ 



(5.5) 



where prime denotes the non- harmonic part of the currents. Thus, only the 
non-harmonic parts of the currents can be bosonized in terms of a single 
valued scalar field. To bosonize their harmonic parts one would have to 
allow for a scalar field 7 with windings as (j) below (3.2). On the infinite 
plane the harmonic part is not present and we may leave out the primes in 
(5.5). If we further assume space time to be flat we recover the wellknown 
bosonization rules in [41] . What we have shown then, is that for the gauged 
model on curved space time the bosonization rules are just the flat ones 
properly covariantized and with the omission of the zero-modes. 

Since (5.4a) holds for any gauge field the current correlators in the 
Thirring model are correctly reproduced by the bosonization rules (5.4,5). 
To see that more clearly we calculate the two-point functions of the auxil- 
iary fields in the Thirring model (5.2,3). For that we differentiate (5.2) {(p 
is treated as external field) with respect to the sources and find 



where Go is the free massless Greenfunction (3.6a,3.8) in curved space-time 
and the integrations are over the variables u and v with the invariant mea- 
sure on the curved torus. Here (. . .)t are vacuum expectation values of the 
Thirring model (5.3a). Alternatively we can calculate these expectation val- 
ues from (4.3) and (4.5), where the fermionic integration has been performed 
and find 




(5.6) 




(5.7) 



Comparing the results (5.6) and (5.7) we see at once that 
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J {Go{x,u)j^^{u)Go{y.v)j'fM)T = 

{Goix,u)j^.^{u)Go{y,v)jl,,{v))r = -^{m'^(p{x)(p{y) - Go{x,y)). 

(5.8) 

The correlators (5.8) express the gauge invariance and the axial anomaly 
Os /i) ~ —fn-yA{p in the gauged Thirring model. They can be correctly re- 
produced with the bosonisation rules (5.5). 

6. Chiral condensate 

We have seen that the partition function of the gauged model does not de- 
pend on the chiral twists p. As we pointed out this property is very much 
linked to gauge invariance. The same should then hold in the topological 
non-trivial sectors. For that reason we shall set /5 = in this section. This as- 
sumption will simplify the manipulations below considerably. For example, 
N^=M^ in (3.34). 

Recalling that in (3.34) anti-commutes with 75 one sees at once that 
only configuration supporting one fermionic zero-mode with positive chiral- 
ity contribute to the chiral condensate 

(V^tp^V) = . r u- r ^ / V{...)ZFMU=n=oe-'-, (6.1) 

where 77+ = P+?7. Earlier we have seen that these are the gauge fields with 
fiux <? = 27r or instanton number k=l. Thus (6.1) reads 

(V'tp^V) = -^/f / 2^(...)V'J+(^)V'o+(x)exp(...)e-^-M, (6.2) 

where exp(. . .) is the last exponential factor in (3.34). Here we have used 
that the chiral twists vanish, such that xo-i- = V'o-i- (see (3.25)). First we 
integrate over the toron field t. The t dependence enters only through the 
zero mode and more specifically ?/'o+ in (3.29) and (3.25) with p — 1. Using 
the series representation for the theta functions one finds 
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J dh^l^{x)iJo+ix) = l. (6.3) 

Note that the result does not depend on the chemical potential similarly to 
our calculation of the partition function. 

To continue we observe that the term J y^G in exp(. . .) vanishes because 
of our conditions (3.36) and (3.38) on the fields ip and 0. Also note, that 
the fermionic Green function does not enter in the expression for the chiral 
condensate. It follows that the fermionic functional (3.33) in the trivial 
sector and (3.34) in the one-instanton sector are the same, up to the factors 
in the fist lines. From (6.3) and (4.3) we see that the toron integral of the 
first line in (3.34) is \ri\^\J tq/V exp(27r J y^x/V^) times the toron integral 
over the factor in (3.33). Also, since 

SB[k=l] = SB[k = 0] + ^ 

the functional integral and normalizing partition function in (6.2) are the 
same, up to these factor and the field-dependent factors which relate the 
hatted and unhatted zero- modes in (3.29). Finally note that the A integrals 
in (6.2) and in the normalizing partition function cancel so that we end up 
with the following formula for the condensate 



4>f 

(6.4a) 



The expectation value is evaluated with 

/-I 2 2 

^\-^{A^ - -A)^ - ^(/>A(/. - ^</.Av. . (6.46) 

A formal calculation of the resulting gaussian integrals yield 



V 



exp [ K[x, x)\ exp [ ^ Gq[x, x)\ , 



(6.5) 



where 



Generalized Gauged Thirring Model on Curved Space-Times 39 
^(^'2/) = (^1 ^2 ,'^^2^ 1^^ " :^{Goix,y)-Gm^{x,y)) (6.6) 

and Gm,Go are the massive and massless Greenf unctions. 

Here we encounter ultraviolett divergences since Go{x,y) is logarithmi- 
cally divergent when x tends to y. To extract a finite answer we need to 
renormalize the operator exp(Q;0). This wave function renormalization is 
equivalent to the renormalization of the fermion field in the Thirring model 
and thus is very much expected [33,41]. In order to do that we first de- 
termine the short distance behaviour of the massless Greenfunction (3.6b). 
Using the idendity 



1 + §1 



2 ^ L 



(0,r)f = |e-(^Wei(^^^,r)|^ 
and the small z expansion 

we see that Go possesses the expected logarithmic short distance singularity 

Goix, y) = -^ log - ^ log (47rVo|77(r)|^) + 0(0- (6.7) 

From the relation (3.8) between the full and hatted Greenfunction and the 
definition of x in (3.4c) it follows that Go has the short distance expansion 

Go{x, y) ^ Go(^, y) + 2x(^) -IJV^X + O(0 (6.8) 



To continue we need to regularize the composite operator exp(Q;0) ap- 
pearing in (6.4a). The normal ordering prescription 
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works well on the whole plane [33,41]. On the curved torus we must be 
more careful when renormalizing this operator. The required wave function 
renormalization is not unique but it is very much restricted by the following 
requirements: First we take as reference system (the denominator in 6.9) one 
with a minimal number of dynamical degrees of freedom since we do not 
want to loose information by our regularization. Second, the renormalized 
operator should have a well-defined infinite volume limit and its expectation 
values should cluster. Finally, the regularization should respect general co- 
variance. These requirements then force us to take as reference system the 
infinite plane with metric g^^. The flat metric is not permitted since it 
leads to a ill-defined expression for (exp(Q;0)). With these choice the nor- 
mal ordering in (6.9) is equivalent to replacing the massless Greenfunction 
in (6.5) by 

Gl^\x,y) := Go(x,y) + ^ log [li's^x.y)]. (6.10) 

Here s{x,y) denotes the geodesic distance between x and y. The occurence 
of the arbitrary mass scale n comes from the ambiguities in the required 
ultraviolett regularization. On the flat torus Gq''^ has now the flnite coinci- 
dence limit 

Gr(.,.) = -^log(l^!l5lf^). (6.11) 

To determine the chiral condensate we also need to determine K{x, y) on the 
diagonal. In a flrst step we shall obtain it for the flat torus. Its cr-dependence 
is then determined in a second step. For a = and t = zto the Greenfunction 
K has been computed in [3]. The generalization to arbitrary r is found to 
be 

miKix, x) = — coth ( , " ) H — 

' 2™,Lr„ |rp ^^^y ^^^^^^ 

+ ^ ( - log I' + F{L, t) - H{L, t)) , 

where we introduced the dimensionless constant a = Lm^\T\/2'n and the 
functions 
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F(L,r) = y[l- , ^ 



1 



1 



+ 



1 



n>0 

We used the abbreviations 



g— 27ri2+(n) g27ri2_(n) 



(6.126) 



= (nri ± %tq\/v?- + a^) . 



(6.12c) 



Substituting (6.12) and (6.11) into (6.5) with cr = we obtain the following 
exact formula for the chiral condensate on the torus with flat metric g/^^: 



-^[F{L,t)-H{L,t))_., 



' coth- 



(6.13) 



where we used that on the flat torus x = and V = V. Futhermore, we 
identified /i with the natural mass scale of the theory. 

To study the finite temperature behaviour of the chiral condensate we 
must assume that t = if3/L and then (3 = 1/T is just the inverse temper- 
ature. Furthermore we perform the thermodynamic limit L — > oo. Then 
coth(. . .) — > 1, iif — > and the expression for the chiral condensate simpli- 
fies to 



(V.tp^V^), = -T(^)^exp 



27r 



(6.14) 



27r + 5r| " 

Let us now investigate the low and high temperature limits in turn. To study 

the low temperature limit we use that 

a 1 

F{I3) ^ ^ + 2 2a a ^ oo, 

where 7 = 0.57721 ... is the Euler number. Inserting this expansion into 
(6.14) yields 

27r 



(V;tp+V^) = -^2^'/(2-+5i) exp ( 



2TT + gi 



7 



for T ^ 0. 



(6.15) 
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For temperatures large compared to the induced photon mass F vanishes. 
Thus we obtain the high temperature behaviour 

(V^tp^V>T = -T(^)^exp(-^T) for T ^ oo (6.16) 

It is instructive to discuss the various hmiting cases. For all gi = 0, i.e. the 
Schwinger model limit, the exact result (6.14) simplifies to 

EtLe^ for T ^ 

47r ^ (6 17) 



where now = c^/tt is the induced photon mass in the Schwinger model. 
This formula for the temperature dependence of the chiral condensate in 
QED2 agrees with the earlier results in [3] . 

Next we wish to investigate how the selfinteraction of the fermions affect 
the breaking. For large coupling g2 and fixed termperature the exponent in 
(6.14) vanishes so that 

{il;^P+tl;)T . ^ for T fixed, g2 00. (6.18) 

V 27r + gl 

Hence, for very large current-current coupling the chiral condensate van- 
ishes. Or in other words, the electromagnetic interaction which is responsible 
for the chiral condensate, is shielded by the pseudoscalar-fermion interac- 
tion. 

For intermediate temperature and coupling g2 we must retreat to numerical 

evaluations of the sums defining the chiral condensate in (6.14). The results 

of the numerical calculations are depicted in Figure 1. 

The study of the influence of the gravitational fleld is complicated by the 

presence of the massive Green function Gm^ in (6.5,6). This Green function 

is known only for very particular curved spaces. Fortunately we only need 

the coincidence limit for which we can use its short distance expansion [42] . 

For simplicity we assume infinite volume and zero temperature. Then [43] 

-,00 
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for small geodesic distances s = s{x,y). Here Hq^^ denotes the Hankel func- 
tion of the second kind and order zero. In particular 

(^)^-[log^+7] for ^^0. (6.20) 



Inserting that into (6.19) we find with Gq — Gq from (6.7) the following 
short distance expansion 

1 r, „,27rWT)|2 



G„(x,,)-G„(x,,)~--[l„g(^mj)-^] 



OO ^ ^ (6.21) 



We have used that ao{x) — 1 and s ~ e"'^^^ s, where s is the geodesic distance 
on the flat spacetime with hatted metric, P = g^i,{x — y)^{x — yY. Finally, 
substituting (6.21) into (6.5) we end up with 

(V^tp^V')a = (V'^J'+V').=o -exp [- i(^)'f;a,(x)^^^^]. (6.22) 

The Seeley-deWitt coefficients aj have been computed up to j = 5 [44] . They 
are of order j in the curvature and its derivatives. The first two are 

ao(x) = 1 and ai{x) — -TZ. 

For TZ « and slowly varying TZ we conclude that the chiral condensate 
decreases with increasing curvature as 

(V^tp^V)~exp[-^7^]. (6.23) 

If wc compare this with the temperature dependence (6.16) we are lead to 
define a curvature induced effective temperature 

For this formal identification of curvature with temperature no horizon is 
needed as in black hole physics where the temperature is related to the 
surface gravity at the horizon. Note that contrary to the temperature the 
curvature may become negative. Then the condensate is amplified. 
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Finally we consider the chiral two point function for non-coinciding 
points. The gauge invariant form reads 

It is related to a bound state between a static external charge and a dy- 
namical fermion [45]. 

The integration over the various fields is similar to (6.3-12). The result takes 
a simple form in the infinite volume and zero temperature limit: 



S+{x,y) = S+{xf^S+{yf^ exp [^{^f {K{x,y) + K{y,x)) 
■ exp [i^^^ + Y^Go{x,y) - ^{Go{x,x) + Go{y,y)) 



(6.25) 



where S+{x) = S+{x,x) — {ip\x)ilj{x)) denotes the chiral condensate. 
Again the massless propagator must be regularized. We do this using the 
prescription (6.10). Then 

, exp [1(^)2 (G',„^,(:r.:^) + G'„.,(:^.:r))" 
S+{x,y) = S+{x)-^S+iy)-^ ^ 



g2 g2 

2^mf^ {g{x)^ g{y)^ s) 

(6.26) 

Note that the coupling strength gi to the longitudinal current enters the 
scaling exponent. On flat space Gm reduces to ■^Ko{ms) wich decays ex- 
ponentially for large separations. Hence we find 

S+{x,y) ^±Ml^ (6.27) 

V27r§{mjs) 

for large separations of x and y. We have used that the chiral condensate 
S+{x) in (6.15) is constant, due to translational invariance. For gi = this 
simplifies to the Schwinger model result [3] 



S+{x.y) - \ rt\ I ^- (6.28) 



2 2'n^J\x - y\ 
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Unlike the correlators of fields which in the bosonized version are local in the 
massive boson field, this two-point function does not decay exponentially. 
However the long range correlations are suppressed by the coupling to the 
longitudinal current. 

7. Thermodynamics of the ungauged Model 

In this chapter we derive the grand canonical potential, equation of state 
and ground state energy for = 0. For the ungauged model there is no 
Gauss constraint and the charge of the vacuum need not vanish. Indeed, for 
Af^ = the partition function depends on the chemical potential and on the 
fermionic boundary conditions. Technically this is due to the absence of the 
toron integration which for the gauged model wiped out any dependence on 
H, a and (3. 

The partition function of the ungauged model is given by 



where Zp is the fermionic generating functional (3.33) and Sb the bosonic 
action (3.35). The integration over the harmonic fields is Gaussian and yields 




(7.1) 




(7.2) 



— oo 



where 




(7.3a) 



is the thetafunction with characteristics 




and 




and covariance 




(7.3c) 
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The remaining functional integrals in (7.1) are performed as in chapter 4. 
To obtain the partition function of the Thirring model in the limit gi = g 

we divide Z by the corresponding partition function Ao of the free bosons 
(see 5.3b)). Using (4.4b) and (3.33) we obtain 



'2n + gi ^[u 



w 



(^)e(^/24-+ff3)5^ (7.4) 



J\fo \viT)\^'\l 27r + gi 

In the Thirring model limit g2 = go and the square-root in this formula 
disappears. 



7.1 Zero-temperature limit 

To investigate the thermodynamics of the model we assume spacetime to 
be flat and that T = i(3/L. Then 

« = -ilog^ (7.5) 

is the grand canonical potential. Let us now investigate the low temperature 
limit of Q. For // = this yields the ground state energy. 
To study this limit we observe that for T = i(3/L the covariance matrix A in 
(7.3c) simplifies to 



iirA = — ^ 
ij 



47r 27r + gl\ -47r - g^ gl 



(7.6a) 



and has eigenvalues 



._l2j^ and A. = -^^^ (7.66) 
L 27r L 2T^ + gl ^ ' 

with corresponding eigenvectors 

vi = (-1, 1) and V2 = (1, 1). (7.6c) 

Also the rj tensor (see 3.17b) and /xq (see 3.11) in (7.3b) simplify to 

. ( (3/L\ , .(3 
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Now we can determine the low temperature limit of the grand potential 
from (7.4) (with Sl = 0) and (7.6). For that we note that the saddle point 
approximation to the Gaussian sum (7.3a) defining the theta-function be- 
comes exact when P —>■ oo. Also, using that 

log|ry(r)p^-|^ for (3 ^ oo (7.7) 

we end up with 

\ Tl" 471 TT uL.2 

«(/?-oc) = ---5^-(ft + |^) 



+ — mm 
2L nez2 



for the zero-temperature grand potential of the ungauged model. The chemi- 
cal potential and chiral twist enter only through the combination (3i+ij,L/27T. 
Up to the second term the potential is invariant under 

Oil — >ai + l and Pi + ^ — > (3i + ^ + 1 + go/^n. (7.9) 

ZTT ZTT 

Let us now discuss the potential in the various limiting cases. 
First assume that there is no chiral twist, Pi — 0, and that the chemical 
potential vanishes. Then f2{P ^ oo) coincides with the groundstate energy. 
The minimum in (7.8) is attained for n\=n2 = [\-\-OLi\ and we find 

Only for anti-periodic boundary conditions, that is for ai = 0, does this 
Casimir energy coincide with the corresponding result for free fermions. For 
^0 — ^Ti" the Casimir force is always attractive whereas for g^ < Att it can 
be attractive or repulsive, depending on the value of cki. The result (7.10) 
is in agreement with the literature [11]. For example, it coincides with De 
Vegas and Destri's result if we make the identification a;^,^, = 27rQ!i and 
1/Pj,j, = 1 + c/2/27r in formula (42) of that paper. 
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For small twists and chemical potential the minimum is assumed for = 
and the potential simphfies to 

and does not dependent on the chemical potential. 

For vanishing go, that is for free fermions, the minimum of (7.8) is attained 
for 

ni = [- + ai - Pi - —\ and n2 = [- + ai + Pi + —\, 
z ztt z ztt 

where [x] denotes the biggest integer which is smaller or equal to x. This 
then leads to the following zero temperature potential 

For = /3i = this reduces to the Casimir energy for free fermions with 
left- right symmetric twists and agrees with the results in [46]. 
Note, however, that for /3i /O we disagree with [11]. The difference is due to 
the second term on the right in (7.8). Let us give two arguments in favour 
of our result: 

The discrepancy arises from the prefactor appearing in the fermionic deter- 
minant (3.19). As discussed earlier this prefactor implies the breakdown of 
holomorphic factorization, a property which has been presupposed in [11]. 
In appendix C we show that our results can be reproduced by starting with 
massive fermions and taking the limit m — > 0. 

The second argument goes as follows: Suppose that j3^—ai = ^. Then (7.12) 
simplifies to 
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For massless fermions the Fermi energy is just /j, and at T = all electron 
states with energies less then fj, and all positron states with energies less then 
—fj. are filled. The other states are empty. Since dQ/d^i is the expectation 
value of the electric charge in the presence of // we see that it must jump 
if crosses an eigenvalue of the first quantized Dirac Hamiltonian h. For 
vanishing twists the eigenvalues of h are just En = {n — \ )'k/L. Indeed, from 
(7.13) one finds that the electric charge 

= ^ = ^^'^ ^""^ ^n<f^<En+l (7.14) 

jumps at these values for fx. Further observe, that in the thermodynamic 
limit L — > cxD the density 

(7.15) 

reduces for 5^0 = to the standard result for free electrons. 



7.2 Equation of state 

We wish to derive the equation of state for finite T in the infinite volume 
limit L ^ 00. This may be achieved by interchanging the roles played by L 
and /3. More precisely, using that 



e 



(yl) = Vdet(iFi) e^^^"^-" © ^" {iA-^) (7.16) 



we find in analogy with the low temperature limit that for L 
pressure is given by 



00 the 



/?p = lim — log 



TT 27r27r+^ 2 
"I ~ Po 



L-ooL ^A/'o 6/3 /? 27r 



TT 

— — mm 
2/3 nez2 



{m +n2 + 2/3o}^ 



27r + (7g 
27r 

21: r ^ „./3/ii2 

+ 9 n2 - ni + 2q;o + 2i— \ 

-•^^ 27r 



27r + ^^ 



(7.17) 
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Here the minimum of the real part has to be taken. Again the minimization 
arises from the saddle point approximation to the theta function which 

becomes exact when L — > cxd. For small twists the miminum is assumed for 
rij = and then 

TT 27r 27r , /3//x2 

becomes independent on the chiral twist Pq. As we have interchanged the 
roles of the temporal and spatial twists this is consistent with the earlier 
result that for small twists f2 is independent of In particular, for ao = 0, 
we are lead to the following equation of state 

TT //^ 27r 

p(/3,;.,a„ = 0) = j^ + ^^, (7.19) 

which for small f3o relates the pressure to the chemical potential and tem- 
perature. This result is consistent with the renormalization of the electric 
charge which is conjugate to the chemical potential. It shows in particular 
that the thermodynamic behaviour of the Thirring model is not just the one 
of free fermions as has been claimed in [12]. Indeed, the zero point pressure 
is multiplied by a factor 2n/{2n + Qq). This modification arises from the 
coupling of the current to the harmonic fields. It can not be seen if only the 
local part of the auxiliary field is considered, which is the case if one quan- 
tizes the model on the infinite Euclidean space. Furthermore, we see that 
the 'pressure' p is real only for ao = 0. This phenomenon occurs also in the 
Hamiltonian formalism [47]. However, finite temperature physics dictates 
anti-periodic boundary conditions, i.e ao — 0, and then p becomes real. 

8. Conformal structure of the ungauged model 

When we discussed the properties of the classical model (1.1) in chapter 2 
we have noticed that for = it reduces to a conformal field theory on fiat 
Minkowski spacetime. We have found the results listed at the end of section 
2. 
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We determine the quantum corrections to these classical results. As in the 
previous chapters we do that within the Euclidean functional approach. 
Thus we start from first principles and need not postulate the emerging 
Kac-Moody and Virasoro algebras in advance [8,16]. When comparing the 
classical with the quantum results one should keep in mind that roles of i/j^ 
and ipl are interchanged when one switches from Minkowski to Euclidean 
spacetime. For further changes the reader is referred to appendix A. 

In what follows it is convenient to exploit the holomorphic structure of 
the model. On the torus with flat metric g^jy the Cauchy-Riemann equations 
read 

(r//a. - za^)/ = 0. (8.1) 

Then one chooses coordinates ^ and the corresponding complex 

coordinates x = x'^ + ix'^ such that (8.1) takes the standard form. More 
explicitly we chose 

X = ifx° + ix^ so that = -^{d^o — Td^i). (8.2) 

2to 

In this section x and x always denote the complex coordinates belonging 
to a;'*. In these coordinates the free Dirac operator and the corresponding 
Greensfunction are simple 

where ^ = x—y. The chiral components of the energy momentum tensor and 
current are then given by 

T.. = |(tT°° + T-) = |^T'^- and j, = l{rf-f). (8.4) 

Using that the energy momentum tensor is conserved and traceless and that 
the vector and axial-vector currents are conserved it is easy to check that 
these chiral components only depend on x and not on x. 
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8.1 Virasoro and Kac-Moody algebras 

First we determine the central charge from the short distance expansion of 
the Txx correlators. As in the classical theory (see (2.13)) the symmetric 
energy momentum tensor measures the change of the effective action F = 
— log Z under arbitrary variations of the metric. For the torus there are 
two independent contributions. One being due to variations of the modular 
parameter r and its conjugate f which depend implicitly on the metric. 
The other is due to the variations of terms which depend explicitly on the 
metric. Since the chiral component Txx is gotten by contracting T^'^ with 
dQuv/df it follows that 

It is always understood when doing metric variations, that we take the flat 
spacetime limit afterwards. The f variation is constant and may be skipped 
in the short distance expansion. 

Taking several metric variations of the curvature dependent part of logZ 
with Z from (7.4), (5.3b) and (4.9) we find the following short distance 
expansions for the three point correlation function 

.rj. rr. rr. . _ 3 + 2477^ 1 

Comparing with the general expression [16] we read off the central charge 
and the conformal weight of the energy momentum tensor 

c = 3 + 2AglTT and Ht,^ = 2 . (8.7) 

The first contribution is that of three free fields. The (73— dependent term 
already appeared in the classical analysis and is related to the coupling to 
the background curvature. It is well known from the minimal conformal 
series. Note that the couplings gi and ^2 do not affect the central charge. 
In particular, if we subtract the central charge of the auxiliary fields and 
set = then the value of c is the same as the one for the Thirring model, 
namely c = 1 [16]. 
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Next we determine the Kac-Moody algebra of the U{1) currents. To de- 
rive the correlation functions with current insertions we couple the fermions 
to a gauge field, that is consider the 'gauged' model without Maxwell term. 
For example, 

Using (4.6) on flat spacetime and without Maxwell term, together with 

Ott 1 

d^ct> = v;Al, where A^^ ^ - —t^ - V ^-W'A, (8.9) 

is the transversal part of one obtains the following short distance ex- 
pansion 

We read off the value k of the central extension in the t/(l)-Kac- Moody- 
algebra to be 



= ^_ , ^3 . (8.11) 



27r 
2TT + g. 

Finally we need to determine the conformal weight of the current. From 



we obtain hj = 1. To summarize, the symmetry algebra is the semidirect 
product of a Virasoro algebra with central charge (8.7) and a U{1) Kac- 
Moody current algebra with central extension (8.11). 



8.2 Conformal weights 

To unravel the possible representations of the Virasoro algebra realized in 
the model we must determine the conformal weights of the fundamental 
fields. The short distance expansions of the fermionic two-point function 
with Tgg follows from the metric variation of the Greensfunction 
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where 



« = 7^-77^)- (8-136) 



and Sij is the fermionic Green function in the external gravitational field 
and harmonic gauge field but with and A set to zero. More precisely, 

However, since Z ~ exp[F(7?.^)], its metric variation vanishes after the flat 
spacetime limit has been taken. We refer to appendix B for the variation of 
Sij and G{x,y). Collecting the most singular terms, we arrive at 



I {J- -J- 

{z — x){z — y) ^z — x z—y' 

2aG{x,y) 



(8.14) 



e 



x — y {z — x)'^ (z — y)^ 2tt z — x z — y 
Using that 

2n x — y 

we find that the 2-point function varies under a infinitesimal conformal 
transformation, paramatrized by f{z), as 

^ dzf(z){Mx) ^l(y) T,,) = [f{x)d, + f{y)dy 

+ ^(1 + ^) [/'(^) + f'iy)] - ^ [/'(^) - f'iy)] } m^miv))- 

(8.16) 

Note that the exponential factor has been absorbed to recover the cor- 
relation function {ilj„{x)iljl{y)). The short distance expansion with T^g is 
calculated similarly. Then one reads off the conformal weights 





1 


= 2 + 1 




= (h^o) 











c-9l i/j_o_ , _2 o 



(8.17) 



1 27rg^ ig^gs 



2 ^"y2 
9i - 



IGtt 2tt + gl 
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Thus we have reproduced the classical results supplemented by additional 
gi and dependent quantum corrections. In the Thirring model limit ^3 = 
and Qi =92=9, these terms add up to give the known anomalous dimension 
appearing in the Thirring model [16]. The last classical term is a peculiar 
feature of the solution. For the conformal weight to be real we are obliged 
to choose ^3 imaginary. 

Let us now turn to the auxiliary fields. It is straightforward to compute 
the correlators 

(Aa; Tzz) ~ -J— 6^3 



Att {x — zY 

(K A. T„) ~ - (8.18) 

{4>x 4>y Tzz) ~ 



IGtt {x — z){y — z) 

We see that the classical results are unchanged, that is for ^3 7^ the scalar 
field A is not primary and for ^3 = we find the conformal weights h\ = h^ = 
0. 

Finally we turn to vertex operators or exponentials of the auxiliary fields. 
In contrast to A and (j) those are well defined even on the extended plane. 
Recalling the regularization prescription (6.10) we find 



1 1 r CKl Q!2 "1^ 



167r 27T+g^ 



+ 

z-x z-y\ (8.19) 



(8.20) 



^. gai<^(a;) . . ga2<^(j/) 

and hence 

- / f{z){: e«^^(-) : : e"^^(^) : T,,) ~ \f{x)d, + f{y)dy 
* Jc '- 

" 8(27rV^|) ^^'-^^^^^ + "^'^f'^y^)] ^""^^'^ ■ ■ ^"'^^'^ '-^ ■ 
From this we read off the conformal weights of the vertex operators 

hi = hi = ' . (8.21) 

8(27r + ^2) ^ J 

Note that ct^ must be imaginary to get a positive weight. A similar analysis 
for the A-field yields 
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- / fiz){: e-^^(-) : : e"^^^^) : T,,) ~ \f{x)d, + f{y)dy 
* Jc '- 



and hence 



(8.22) 



(8.23) 



Here both ai and ^3 must be imaginary for the weights to be positive. Note 
that contrary to A the fields : e"^'^^) : remain primary when the XTZ coupHng 
is switched on. This coupling results only in a shift of the conformal weights. 



8.3 C7(l)-charges 

To see how the left and right Kac Moody currents act on the fermionic 
fields we notice that after the integration over the auxiliary fields the A- 
dependence of the fermionic Greenfunction factorizes as 



where 



g{x) = -ia{x) + 75/3(^(a;), P 



2tt 



27r + ^i 

Also, using that on fiat spacetime 

0(a;) = -iJ d^G{x, z)A^ + i J d^G{x, z)A 
a{x) = J d^G{x, z)A^ + j d^G{x, z)A^, 

one ends up with 

■An + gi 1 (7? 1 



(8.24) 



1 

1 
47ri 





2tt + gl z-x 27r + 5f| z-y 

,2 



9l 



2-K + gl z-x ' 2-K + g^ z-y 



1 ^Aii + gi 1 



{M^)iPi{y)^) 

(8.25) 



and thus obtains the following the ^7(1) charges 



1 



27r 



1 



27r 



27r + 5r. 



(8.26) 
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We have used the convention where the electric charge q+q is unity. In the 
Thirring model limit we can compare (8.26) with the results obtained in 
[16]. For that we need to rescale the currents such that the central extension 
(8.11) of the Kac-Moody algebra becomes unity 

j,^Vl + g^j27rj,. (8.27) 

Now it is easy to see that we agree with [16] if we make the identification 

g^^ = ^ ^ = . (8.28) 

Let us summarize our results. The classical conformal and axial transfor- 
mations of all fields besides </> and A are deformed. The longitudinal part 
of the current-current interaction in (1.1) changes the conformal weights of 
the fermion field only. The transversal part affects all weigths and U{1)- 
charges. The background charge changes the conformal weight of the vertex 
operators belonging to the scalar field. 

Of course, the same structure is found in the other chiral sector. 
8.4 Finite size effects 

When quantizing a conformal field theory on a spacetime with finite volume 
one introduces a length scale. The presence of this length scale in turn breaks 
the conformal invariance and gives rise to finite size effects. It has been 
conjectured [13] that the finite size effects are proportional to the central 
charge. For example when one stretches space time, — > ax", then the 
change of the effective action is proportional to c: 

Faa^-r^^—loga-x, (8.29) 

where x is the Euler number of the euclidean space time. In [48] this conjec- 
ture has been proven for a class of conformal field theories on spaces with 
boundaries. The only important assumption has been that the regurlariza- 
tion respects general covariance. In this subsection we shall show that the 
conjecture does not hold for the model (1.1) on Riemanniann surfaces. 
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Unfortunately, the only global conformal transformations on the torus 
are translations which do not give rise to finite size effects. Also, the Euler 
number vanishes and according to (8.29) the finite size effects are insensitve 
to the value of c. For that reason we quantize the ungauged model (1.1) on 
the sphere where the global conformal group is the Moebius group. 

An effective method to compute finite size effects has been developped in 
[48] . It is based on the following observation: Any conformal transformation 
z w{z) is a composition of a diffeomorphism (defined by the same w) and 
a compensating Weyl transformation g^jy —>■ e^^ g^i, with 

,2. ^dwiz^dMz^^ . = + (8.30) 
dz dz 

Therefore, chosing a diffeomorphism invariant regularization one has 

= SFoiff = SFconf - ^Fweyl- (8.31) 

Now we apply the techniques of the previous sections to derive the change 
SFweyi of the effective action on the sphere under Weyl transformations. 
This change is given by the trace anomaly. 
The change of the effective action under Weyl rescalings is 

Sr„,,. = - log S^i^4.)ietmeM-SslA = o,s]) 

JV{X(t))det{ip)exp{-SB[A = 0,g]) 

where Sb is the bosonic action (3.35) with vanishing gauge field. Also, since 
on the sphere there are no harmonic vector fields the term ~ h'^ in Sb is 
not present. Thus the calculation on the sphere is actually simpler as on 
the torus (see 7.1) since there is no integration over the harmonic fields. 
As on the torus we must impose the conditions (3.38) in order to eliminate 
the additional degrees of freedom we introduced in the derivative coupling 
representation. Thus we obtain 
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Here we used that (3.30a) in the trivial sector still holds on the sphere. Also 
we used the scaling law (4.4b). Sl is the Liouville action (3.30b) in which 

we can not put TZ to zero, since 



/ 



= Stt = 47rx (8.34) 



for any curvature and thus in particular for As for the fermions (see 3.26) 
one introduces the 1-parametric family of Laplacians 



= e-^^^A (8.35) 

interpolating between A and A. The r derivative of the corresponding de- 
terminant is given by the trace anomaly [32,48]. The explicit calculation 
yields 

logdet'^ = 2 j^dr j - — < - (8.36) 

Again g'^ is the determinant and a[ = ^TV the relevant Seeley-deWitt 
coefficient of the deformed metric g'lij^ = e^'^'^ g^v is the projection onto 
the zero- mode of A"^. Using that the normalized zeromode is constant and 
~ 1/ \/V^, one finds 

logdet'^=log^ + -^5z,. (8.37) 

A V J-^TT 

The ~ log V term cancels against the same term in (8.33) and we end up 
with 

5r = 4 / VgnUn -^^) - ^ [ ^gna+ ^ [ ^gaha. (8.38) 



4 J ^ A ' V ' 2Att J ^ 247r _ 

Now we can see why the finite size conjecture generally fails to be true, 
although it holds for theories without background charge on domains with 
boundaries [48]. Take the simple case of a dilatation w{z) =az. Then, the 
conformal angle is a constant cr = loga and {TZ — 87r/F) = 0. Thus the first 
term in (8.38) vanishes and the finite size effect does not depend on g'^. It 
is given by 

6r = -^loga I V^7t=-loga (8.39) 
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and does not agrees with (8.29) since c in (8.7) depends on g^. Thus we 
have disproved the conjecture. On other Riemannian surfaces one would 
find the same result: the effective action scales as in (8.7) where c is the 
central charge of the model without background charge. It is evident that 
the finite size scaling comes from the middle term ~ logo J \/gR in (8.38). 

It is interesting to compare the finite size scaling on Riemannian surfaces 
with the one on domains with boundaries. In the presence of boundaries 
(8.36) is modified to 



where the second integral is over the boundary of spacetime and g^^ the 
induced metric on this boundary. On a domain we can always put R to 
zero and the middle term in (8.40) does not contribute to the scaling. The 
scaling comes from the surface term in (8.40). Diffeomorphism invariance 
implies that the bulk term determines the surface term (up to diffeomor- 
phism invariant surface terms). This is how the central charge, defined by 
the short distance expansion of the Tgz-correlators and thus by the bulk 
term, reemerges in the scaling law (8.7), which is determined by the surface 
term. 




(8.40) 



Acknowledgements 



This work has been supported by the Swiss National Science Foundation. 
We would like to thank K. Gawedzki, C. Kiefer and E. Seller for discussions. 

Appendix A: Conventions 

In this appendix we set up our notation and give a list of useful formulae. 
Let g^i, be the metric of spacetime. The sign convention for the curvature 
tensors is such that 

In 2 dimensions the only independent component is TZowi - In order to couple 
fermions to gravity we must introduce a local Lorentz frame (or tetrad), e^a, 
relating the Lorentz and spacetime indices: 

^aa^/s' = dap 1 ^aa^\ ~ Vab 1 Vab ~ (^Q ~1 ^ (^•^) 

The latin and greek indices are Lorentz and spacetime indices, respectively. 
All physical laws should be general- and Lorentz covariant. If gai3 has Eu- 
clidean signature then r]ab in (A. 2) is changed to dab- 
The 'curved' gamma matrices are related to the flat ones as 

We us the following chiral representation for the flat 7's: 

and in Euclidean spacetime we may choose 
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We may also define 
The relations 

7^75 = 6'';,7' , l%% = -ie\l\ where eab=(^\ {AM) 

are particular to 2 dimensions and play an important role in the body of 
the paper. Note that depending whether one is in Minkowskian or Euclidian 
spacetime the Lorentzindex a is raised with ■q"'^ or S"-^. The curved space 
analogue of (A. 6a) reads 

75 = \r],xul^lM = Yi^nul^lE = 75, (^-7) 

where rj^v = \/\9\^^.v is the antisymmetric tensor (whereas the flat metric has 
Lorentz- indices, the antisymmetric tensor has space-time indices). To im- 
plement local Lorentz invariance one needs to introduce a connection a;^o6- 
For example, in the Lagrangean the Lorentz-covariant derivative acting on 
the spinors read 

D^, = d^ + iuj^, (A.8) 
where the spin connection uj^ is defined by 

(A.9) 

In 2 dimensions this reduces to 

1 1 

Finally we list some useful scaling relations. If the 2-bein scales as — e'^e^ 
then the above introduced quantities scale as 
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Appendix B: Variational formulae 



In the following denotes the spacetime and Lorentz covariant derivative. 
How it acts on spacetime and Lorentz tensors follows from the first formula 
in (A.9). 

Using the definition of the Christoffel symbol and (A. 2) it is straight forward 
to show that 

Sg,j,u = ^e^eya + e^be^a , S^/g = ]^y/gg'^''5g^^ 

Sr^. = Ig^'^iDJgp^+D^Sgf,, - Df,5g^,). 
For some formulae related to the variation of the tetrad let us refer to [49] 

{B.2) 



where t\ = ^(e^«<5e,6 - e^^^e/). 
Then using (A.9) it is easy to see that 

When performing the variation of curvature dependent expressions we have 
used the identities 

g^^''Sn^.=u%, where = g^'^SF^, - g^'^SFi:, 

and J ^Lo'^A^ = J ^{g^f'V^A'' - V«A^}5^a^ . ^^'^^ 
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Depending on the topology of spacetime, the induced curvature 1Z appearing 
in (A. 11) may be different from zero. In this case it is not possible to express 
the conformal angle a in terms of the curvature scalar. Nevertheless, to 
perform variations of d-dependent expressions, the identity 

proves to be useful. 

Taking the variations of the equations 

^AG{x,y) = -5{x-y) and ^iIpS{x,y) = 5''{x - y) (B.G) 

for the scalar and fermionic Greensfunctions we may derive (up to contact 
terms) the following variational formulae 

SS = '-J {2S{x, u)rD''S{u, y) -D^[S{x, u)^%^v''"S{u, y)]) ^dg^., 

{B.9) 

where all arguments and derivatives which are not made explicit in the 
integral refer to the coordinate u over which is integrated. Finally, we need 
the following formula for the variation of the inverse Laplacian 

* (z^) = M*-^ " **^'^-^) " ^ / ^^""^^""^f' ('^■^°) 

where V is the volume of spacetime and / an arbitrary function. To prove 
this identity we note that for / e (KernA)-*- we have 

Varying this equation yields 



A(5^/) = 5/-(5A)l/ 



which may be inverted to give 
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Varying the identity 

If 1 

allows to replace the last term of (B.ll) to obtain the required result (B.IO). 

Appendix C: Canonical approach to the partition 
function 

In this appendix we compute the partition function for massive Dirac 
fermions in the canonical formalism. In the limit m ^ we confirm some 
of the results in sections 3 and 7. For massive fermions one cannot con- 
sistently impose chirally twisted boundary conditions. However, from the 
explicit eigenvalues (3.17) one sees at once that the chiral twist /3i and the 
chemical potential are equivalent. One can easily verify that this equiva- 
lence holds also for massless fermions in the canonical approach and that 
Pi ~ (jlL /2tt. Let us therefore compute the partition function 

for massive Dirac fermions with chemical potential on a cylinder with 
(nonchiral) twisted boundary conditions 

■4}{x + L,t)^ -e-^'^^^'tPix, t). {C.2) 

For massive particles it is more convenient to use the Dirac representation 

/ = a3 7^ = -z(72, 7^ = 7°7i = -(71. (C.3) 

The Dirac field is expanded in terms of the eigenmodes of the first quantized 
Hamiltonian 




(CA) 



as 

W{X, t) = J2 '^n, + hn + V'n,-<, (C.5a) 
n n 
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where the V'n,+ ^ind V'n,- ^ire the positive and negative energy modes, 

Un + m\ (C.56) 



-'n 



^/2uJn{cOn +m)L V 

The momenta and energies are determined by the boundary condition 
(C.2) to be 

27J- 1 

Xn = Y^{n - - - ai) and u>n = \J ^ \\. (C.5c) 

After normal ordering the 'positron' operators with respect to the Fock 
vacuum defined by H we find 

n n n 

where the last c-number term represents the infinite vacuum contribution 
which must be regularized. To do that we employ the zeta function regular- 
ization. That is we define the zeta-function for s>l by the sum 

n 

which in turn defines an analytic function on the whole complex s-plane 
up to a simple pole at s = l. The analytic continuation is constructed by a 
Poisson resummation 



n 



where 



and rh = Lm, fi = Lji. Taking the Mellin transform of (C.7b) we find 

=;^e^-^(^-«^) j dt t^-'e-'f^ j^Ko{fi^/e+^) (C.7c) 
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F diverges at since the Kelvin function Ki{z) ~ 1/z for small z. It follows 
that the n = term in (C.7a) diverges. This divergence is regularized by 
substracting the ground state energy of the infinite volume system. Indeed, 
because of the exponential decay of the Bessel function for large arguments, 
only the n = term contributes for infinite volume. So we find for the 
regularized sum 



Noe we perform the limit m — 0. Only the most singular term in the 
expansion of the Bessel function contributes, hence 



where Sa-biz) is the Lommcl function [50]. In particular for s = —1 this 
function is S = 1/z so that finally 

Inserting this into (C.6) then yields the regularized expression 



: H-nQ := ^{uJn - I^Wnbn + ^{uJn + I^Wndn - — + j^i^^ - 1.^^ + 2^)'^ ■ 
n n 

(C.ll) 

For small the normal ordering is //-independent so that 

(0| : iy - /.g : |0) = + ^(ai - [a, + ^])^ - (0| : if : |0) {C.12) 



is independent of and coincides with the Casimir energy [46] . 

Finally we compute the partition function. Using (C.12) we easily find 
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Z{P) = tr [e-l3■■iH-^^Q)■■] = gl^l-^^] 

oo 



n (1+^ 



(n-i+ai)g/3/x^_ 



n>[|+ai] 
oo 

n (1+'! 

n>[^-ai] 



n>-[^+ai] 
oo 



(n-i+ai)g-^/u 



n>-[i-ai] 



— tti 



(O,r)0 



— cti 



](0, 



T 



(ai3) 

where we have used the product representation of the theta functions in the 
last idendity and that q = e^'^*'^ = e"^'^^/-'^. A non- vanishing chiral twist 
Pi can now we included by shifting the chemical potential. Thus we have 
confirmed the formula (3.19) in our paper. 

Note that for n ^ the zero-temperature limit of the grand potential is 
not equal to the vacuum expectation value of : H — jiQ :. For 7^ all states 
up to the /i-dependent fermi energy are filled. For example, for u}\ < jj. < u}2 
the limit lim/3^00 ^ reduces to the expectation value of : H — fj,Q : in the 
one-electron state. 
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